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Abstract 

Large-scale antenna (LSA) has gained a lot of attention due to its great potential to significantly 
improve system throughput. In most existing works on LSA systems, orthogonal frequency division 
multiplexing (OFDM) is presumed to deal with frequency selectivity of wireless channels. Although 
LSA-OFDM is a natural evolution from multiple-input multiple-output OFDM (MIMO-OFDM), the 
drawbacks of LSA-OFDM are inevitable, especially when used for the uplink. In this paper, we 
investigate single-carrier (SC) modulation for the uplink transmission in LSA systems based on a 
novel waveform recovery theory, where the receiver is designed to recover the transmit waveform while 
the information-bearing symbols can be recovered by directly sampling the recovered waveform. The 
waveform recovery adopts the assumption that the antenna number is infinite and the channels at different 
antennas are independent. In practical environments, however, the antenna number is always finite and 
the channels at different antennas are also correlated when placing hundreds of antennas in a small area. 
Therefore, we will also analyze the impacts of such non-ideal environments. 
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I. Introduction 

Large-scale antenna (LSA) has gained a lot of attention recently [Hl-01. In an LSA system, 
a base station (BS) is equipped with hundreds of antennas. It can be considered as an extension 
of the traditional multiple-input multiple-output (MIMO) system, which has been widely studied 
during the last couple of decades 0. The properties of LSA make it a potential technique for 
future wireless systems. 

Through the employment of a large number of antennas at the BS, the channel vectors for 
different users will be asymptotically orthogonal if the channels corresponding to different users 
and different antennas are independent 0. In this case, the matched filter (MF) becomes the 
optimal detector J6l. The asymptotical orthogonality of channel vectors corresponding to different 
users also allows multiusers to work in the same bandwidth, and thus can improve the spectrum 
efficiency of the network 0. It is shown that the transmit power for each user can be scaled 
down by the antenna number or the square root of the antenna number, depending on whether 
accurate channel parameters are available or not 0. As a result, the transmit powers of users 
can be arbitrarily small when the number of antennas is large enough, and the energy efficiency 
of users can be therefore significantly improved. 

Traditional MIMO systems are usually combined with orthogonal frequency division multi¬ 
plexing (OFDM), where the latter is used to convert frequency selective channels into flat fading 
channels. As a natural evolution of MIMO-OFDM, OFDM has been presumed in most existing 
works in LSA systems 0, flU, [0. Although straightforward, the drawbacks of LSA-OFDM are 
inevitable, especially when used for the uplink. The high peak-to-average power ratio (PAPR) 
of OFDM will reduce the efficiency of the power amplifier at the user terminal, making OFDM 
not suitable for uplink transmission. Meanwhile, the BS will need hundreds of discrete Fourier 
transforms for OFDM demodulation, resulting in a heavy computational burden. On the other 
hand, it is well known that single carrier (SC), due to its low PAPR, has replaced OFDM for 
the uplink of long-term evolution (LTE) even though the latter has been widely advocated for 
the downlink. Following the same philosophy as LTE uplink, SC should be still favorable for 
the uplink in LSA systems to save the needs of high performance power amplifiers at the user 
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terminals. In view of the advantages of SC over OFDM when used for the uplink, LSA-SC is 
expected to be a competitive candidate for uplink transmission in future wireless systems. 

Due to the reasons above, SC has been investigated in LSA systems, especially for the uplink. 
In (91, we have considered SC for the uplink transmission in an LSA system over Rician fading 
channels. In that case, an equalizer is required at the receiver to suppress multiuser interference 
caused by the line-of-sight path. In IITOll . SC is used for precoding with an MF precoding matrix. 
The work in iflOl can be easily extended to uplink according to the duality [0. Traditionally, 
the optimal multi-antenna SC receiver is composed of a set of analog MFs followed by an 
equalizer that is used to suppress the inter-symbol interference (ISI) [fTTI - |[T3l . The need of the 
equalizer significantly increases the complexity of the SC receiver due to the matrix inverse. 
This is especially the case when SC is used for high-speed data transmission. 

In this paper, we propose an LSA-SC receiver where the ISI is suppressed through the LSA 
array. The LSA-SC receiver is only composed of a set of analog MFs and the equalizer is not 
required any more. Hence, the complexity of LSA-SC receiver can be greatly reduced compared 
to the traditional multi-antenna SC receiver ifTTI . Ifl2l . Essentially, the LSA-SC receiver in this 
paper can be viewed as a special case of the traditional multi-antenna SC receiver when the 
antenna number is infinite and the channels at different antennas are independent IfTTI . lfT2l . As 
an alternative, our derivation is based on a novel waveform recovery theory, where the receiver 
is designed to to recover the transmit waveform while the information-bearing symbols can be 
recovered by directly sampling the recovered waveform. Compared to existing works IfTTI . Iff2ll . it 
provides a much simpler way to understand the behavior of SC receiver in large-scale regime. On 
the other hand, the waveform recovery theory adopts the assumption that the antenna number 
is infinite and the channels at different antennas are independent. In practical environments, 
however, the antenna number is always finite and the channels at different antennas will be 
correlated when placing hundreds of antennas in a small area. We will therefore analyze the 
impact of the non-ideal environments in this paper. 

The rest of this paper is organized as follows. Section II introduces the system model. In 
Section III, we investigate the waveform recovery in an LSA-SC system. In Section IV, we will 
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discuss the impact of non-ideal environments on the LSA-SC receiver. Numerical results are 
presented in Section V. Finally, conclusion is drawn in Section VI. 


II. System Model 


In this section, we will first introduce the channel model, and then we will present SC 
modulation for an LSA system. 

A. Channel Model 

For a frequency-selective fading channel, the basedband channel impulse response (CIR) at 
the m-th receive antenna can be given by 



( 1 ) 


where a m [l\ is the complex gain of the Z-th tap on the m-th antenna, 77 is the corresponding 
tap delay, and <5(-) represents the Dirac delta function. In (Q]), we have implicitly assumed that 
the tap delays for different antennas are the same since the extra propagation delays caused 
by the physical size of the antenna array are quite small and thus can be omitted. If the gains 
corresponding to different taps are independent and complex Gaussian, then 


E {a m [l}a* m \p}} = of 6 [l — p\, 


( 2 ) 


where of is the power of the Z-th tap and £[•] denotes the Kronecker delta function. 


Denote c(t) = [ci(Z), • • • ,c M (t)] T to be the CIR vector corresponding to different antennas 
at the BS. From (OQ), we have 



(3) 


where ot[l] = (a^l], ■ ■ ■ , %[!]) T denotes the complex gain vector corresponding to the Z-th tap. 
If assuming the antenna number in an LSA system is infinite and the complex gains at different 
antennas are independent, the mean in ([2 ]) is then equal to the sample mean according to the 
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Fig. 1. SC transmission with multiple receive antennas. 


law of large numbers lfl5Tl - that is 

I™ j7Ct n [l]a\p] = V {am[l}a* m \p}} 

M—>■ oo lVl 

= a?5[l-p\, (4) 

since cq[m]’s with m = 1, 2, • • • , M are independent and identically distributed random variables. 
The last identity indicates that the complex gain vectors corresponding to different taps will be 
asymptotically orthogonal. 

B. SC Modulation 

Fig. [H shows the uplink SC transmission in an LSA system with M receive antennas. Without 
loss of generality, we assume that the user is equipped with only one transmit antenna for 
simplicity even though our result can be directly applied to the multiple antenna case. If the 
transmitted symbols, s[/c]’s, are independent with zero mean and unit variance, then we have 
E{s[£;]s*[n]} = 5[k — n]. Given a transmit filter, g T (t), the transmit waveform is 

OO 

s T (t) = s[k\g T (t — kT), (5) 


where T is the symbol period. 
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Denote r (t) = [r | (t), ■ ■ ■ , r M (t)\ v to be the received signal vector at the BS, then we have 

r(f) = ^2 a [ l ] s T(. t ~ T i)+ v(t), ( 6 ) 

i 

where v(t) = [iq (t). ■ ■ ■ ,%(i)] T is the additive white Gaussian noise (AWGN) with a constant 
spectral density, N 0 . After going through a receive filter, g R (t), at each antenna, the filtered 
receive signal vector, x(t) = [x\ (t), • • • ,XM(t)] r , becomes 


x(t) = r (t) * g R (t) 

= ^ a [l}s(t - n) + z (t), (7) 

i 

where s(t ) is given by 

OO 

s(t) = s i k \ 9 (t~kT), (8) 

k =—oo 


with p(f) = g T (t) *g R (t) and z (t) = [zi(t), ■ ■ ■ , z M (t)] T is the filtered noise vector with z m (t) = 

Vm(t)*g R (t). 

In general, the same root-raised-cosine function with roll-off factor, (3 , is used for both transmit 
and receive filters, and the overall impulse response thus becomes a raised-cosine function which 
can satisfy Nyquist condition liT4ll . that is, g(kT) = S[k] and thus s(kT ) = s[k\. In this case, 
the correlation matrix of the filtered noise can be expressed by 


E{z(t)z n (t + r)} = N 0 g(r)I, 


(9) 


where I denotes the identity matrix. 


III. LSA-SC System 

In this section, we will first introduce the single-tap LSA-SC receiver using waveform recovery, 
then a multi-tap receiver will be addressed. Finally, the relationship with the traditional multi¬ 
antenna SC receiver will be discussed. 
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Fig. 2. LSA-SC receiver with (a) single-tap (b) multi-tap. 


A. Single-Tap Receiver 

As shown in Fig. [2] (a), a single-tap receiver exploits a particular tap from the received signal 
to recover the desired waveform, s(t). To recover the waveform over the p-th tap with delay r p , 
as in Fig. [2] (a), we have 

s P (t) = -^a H [p]x(t + t p ), (10) 

Using the asymptotical orthogonality in ©, the complex gain vectors are asymptotically or¬ 
thogonal if assuming the antenna number is infinite and the channels at different antennas are 
independent. In this case, (fTOl) can be obtained as 

s p (t) = o%s(t) +w p (t), (11) 

which is composed of the desired waveform, s(t), except for being scaled by a factor, a 2 , and 
the additive noise, w p (t), with 

w p {t) = -^c* H [p]z(t + r p ), (12) 

whose correlation function is given by 

E {w p (t)w*(t + r)} = ^jg(r)a 2 p . 


(13) 























From dTTT) . the multipath propagation has been completely removed using the asymptotical 
orthogonality in ©. The desired waveform can be therefore recovered over the tap of interest 
except for a scaling factor. 


Note that the waveform recovery theory only provides a way to understand the SC receiver in 
large-scale regime. The purpose of a communication system is still to recover the information¬ 
bearing symbols. In this sense, a sampler is required for the recovered waveform, as shown in 
Fig. [2] (a). Since g{t) satisfies the Nyquist condition, the symbols can be easily recovered with 
baud-rate sampling. Therefore, the decision variable corresponding to the p-th tap at t = kT is 


'Spik] = cr 2 s[k] + w p (kT). 


(14) 


Accordingly, the signal-to-noise ratio (SNR) corresponding to the p-th tap will be 



SNRp = a 2 


(15) 


B. Multi-Tap Receiver 

In above, the desired waveform is recovered with only the p-th tap. Actually, the SNR can 


be further improved by exploiting all taps through a linear combination of s p (£)’s with p 


0,1, • • • , L, as shown in Fig. [2] (b). 

Let r/ [p] with p = 0,1, ■ • • , L be the coefficients for the linear combination, then the overall 
recovered waveform from all taps can be expressed as 


L 


s(t) = ^2ri\p]s p (t) 


p= o 


L 


L 



( 16 ) 


p =o 


p =0 


Similar to the single-tap case, the decision variable after sampling at t = kT can be obtained as 


L 


L 


s[fc] = J ^r][p}a 2 p s[k ] + ^ rj\p]w p (kT). 


P =o 


p= 0 


(17) 


Accordingly, the SNR after combination will be 
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SNR = 


(18) 


p=0 

Direct calculation from (fl8l) yields that the maximum SNR can be achieved when 


V [Q] = v [l] = ■■■ = r)[L\, 


(19) 


and the maximum SNR, if we assume that the total channel power is unit, will be 



( 20 ) 


for any p. As expected, the SNR of the multi-tap receiver can be improved compared to the 
single-tap receiver. 

C. Relation to Traditional Multi-Antenna Receiver 

Without loss of generality, we assume that rj\p] = 1 for simplicity. Then, by substituting (fTOl) 
into (fl6l) . we can obtain 



( 21 ) 


where © has been used for this equation. The receiver structure for (|2TI) is shown as in Fig. [3] 
From the figure, the LSA-SC receiver is only composed of a set of analog filters and the equalizer 
is not required any more. Hence, the complexity of LSA-SC receiver can be substantially reduced 
compared to the traditional multi-antenna SC receiver ffl2l . 

Essentially, the LSA-SC receiver can be viewed as a special case of the traditional multi¬ 
antenna SC receiver if the antenna number is infinite and the channels at different antennas are 
independent. In that case, the equalizer in traditional multi-antenna SC receiver will reduce to a 
Kronecker impulse, and thus can be omitted in the receiver. Although the resulted receivers are 
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LSA-SC receiver 


Fig. 3. From traditional multi-antenna SC receiver to LSA-SC receiver. The Equalizer reduces to a Kronecker impulse and 
thus can be omitted in the LSA-SC receiver. 


similar, the waveform recovery theory based on the asymptotical orthogonality in LSA systems 
provides a straightforward way to understand the behavior of SC in large-scale regime. This 
is especially the case when the LSA-SC receiver works in non-ideal environment as it will be 
discussed in Section IV. 


IV. Impact of Non-Ideal Environments 

The derivation of the LSA-SC receiver in Section III assumes that the antenna number is 
infinite and the channels at different antennas are independent. In practical systems, however, 
the antenna number is always finite and the channels at different antennas will be also correlated 
when placing hundreds of antennas in a small area. Therefore, we will investigate the impact of 
non-ideal environments on the LSA-SC system in this section. 
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D 



m — 0 1 M — 2 M — 1 


Fig. 4. ULA: antennas are placed along a line with length D. 


A. Spatial Correlation 

Without loss of generality, we consider a typical uniform-linear-array (ULA), where M 
antennas are placed along a line with length D as in Fig. |4j The spatial correlation analysis 
here is also valid for any other form of antenna arrays, and the correlation function will then 
depend on the geometry of the array. In the ULA, the distance between the m-th and the mo-th 
antennas will be d = (m — m 0 )D/(M — 1). Denote the spatial correlation between any two 
antennas with distance, d, to be 

The particular form of /?(•) depends on the surrounded environments. If the power azimuth 
spectrum is a constant from 0 to 27r [fI71 . then 

p(d/X) = J 0 (2nd/X) 1 (23) 

where J 0 (-) is the zeroth-order Bessel function of the first kind. 

In the presence of spatial correlation, we have 

E {]} = o]p ' j ) W-p]- ( 24 ) 

Accordingly, the correlation matrix for the tap vectors, a[l], is given by 


E {a[(]a H [p]} = af6[l-p\ R, 


(25) 
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where R is the spatial correlation matrix, with (m, m 0 )-th entry being 

(m — mn D\ 

[ R ]( m ' m ») = P \ M - 1 ' ~XJ ‘ ^ ^ 

Since the channels at different antennas are correlated, the asymptotical orthogonality in © is 
not valid any more. 

B. Residual ISI 

If the antenna number is finite or the channels at different antennas are correlated, the recovered 
waveform can be obtained, from (I2TI) . as 

s(t) = Jf^2\\ a \p]\\l s(t) + ^J2J2 aH ^ a ^ t ~ Tl + T P^ + ^2 w P^ t ^ (27) 

p p p¥=i p 

where the first term in (|27T) is the desired waveform, and the second term is the inter-wave 
interference (IWI) which denotes the interference to the desired waveform caused by the multi- 
path propagation. If the antenna number is infinite and the complex tap gains corresponding to 
different antennas are independent, the IWI can be completely suppressed using the asymptotical 
orthogonality, as we have discussed in Section III. However, it cannot be completely removed in 
the case of finite antenna number or in the presence of spatial correlation, resulting in residual 
ISI to the decision variable after sampleing. Note that the IWI is in general different from the 
ISI. They are only equivalent when tap delays are integer times of symbol period. 

To analyze the residual ISI, we rewrite the received symbol, from Fig. [3} in a symbol spaced 
form as 


s[k ] = f[0)s[k] + ^ f[n\s[k — n] + w[k\, (28) 

njt 0 

where f[n] is the overall discrete impulse response observed by the receiver, that is, 

/n=4ee <x H \p\ct[t\g(nT-Ti + Tp), 

i p 


(29) 
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and w[k] is given by 

w[k] = ^^Wp(kT). (30) 

p 

In m, the first term is the desired symbol, the second term is the residual ISI, and the third 
term is the additive noise. 

In the ideal case, that is, the antenna number M —> oo and the channels at different antennas 
are independent, the complex gain vectors for different taps will be asymptotically orthogonal 
and thus f[n] = S[n]. As a result, the ISI in (l28l) can be completely suppressed. 

In the non-ideal environments, the asymptotical orthogonality in © is unavailable. In this 
case, f[n] ^ 6[n] and thus ISI arises. The residual ISI can be expressed as 

I[k] = ^2 f[n]s[k - n]. (31) 

n^k 

Direct calculation from the Appendix shows that the average power of residual ISI can be 
expressed as 


where 


j=.s. = E{|/[fc]| 2 } 

= Po YiY a l a 'is 2 ( nT - t i+ t p ). 

n^O p l 


Po = 


tr{R 2 } 
M 2 


(32) 


(33) 


In (l32l) . Pisi is determined by the spatial correlation through P 0 and the power delay profile. To 
get more insights, we consider the following three different cases. 

1) Finite Antenna Number: In this case, we have R = I since the channels at different 
antennas are independent, but P 0 = 1/M rather than zero since the antenna number is finite. As 
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a result, (l32l) reduces to 


p ^i = jj^2J2J2 a2 p^ g2 ^ nT ~ ti + t p^ (34) 

0 p l 

From (l34l) . the power of residual ISI is inverse proportional to the number of antennas. In an 
extreme case, we have Pjsi = 0 if M —> oo, which agrees with our analysis in ideal environments. 

2) Spatial Correaltion: If the channels at different antennas are correlated, tr{R 2 } can be 
rewritten, from the Appendix, as 

tr { R2 }= p2 [ J n _ i 'y) ( M ~ H)- (35) 

m=—(M—1) ' ' 

Therefore, when M — » oo, P 0 can be rewritten as 

Po = I P 2 (jjp) “ \ x \) dx > 0 - ( 36 ) 


It means that in the presence of spatial correlation, the ISI cannot be completely suppressed 
even when the antenna number is large enough. 

In addition, direct calculation of (l36l) for the typical Jakes correlation in (1231) yields that |[T8l 


-Po — 2 • 2P3 
~ 2 • 2P3 




4tt 2 T) 2 \ 

4tt 2 T) 2 \ 

) 
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(37) 


where 2 -P 3 (-;*;■) is the generalized hypergeometric function and ./](•) denotes the first order 
Bessel function of the first kind. The last identity is due to the fact that ,/ (l (•) and ,/](•) are very 
small when D/X is large enough. Since 2 Fs(-', •; •) increases as the decreasing of D/ A, the ISI 
power can be therefore reduced by using longer ULA. 

3) Integer Tap Delays: In this case, the tap delays are integer times of symbol duration, that 
is, Ti = IT. Then, we have g 2 (nT — T/ + r p ) = 8[n — l + p], and substituting this identity into 
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, we can obtain 


Pisi — Po (1 - ^2 


cr, 


(38) 


1=0 


Apparently, P ISI in (l38l) can achieve its maximum when erf’s are all equal if assuming a i — 1- 
In this case, we have 


P 


ISI,max 


L + l 


Po ~ Po, 


(39) 


where the second approximation holds on when L is large. It means that the maximum of the 
residual ISI power is irrelevant to the power delay profile but entirely determined by the spatial 
correlation 


V. Simulation Results 

In this section, numerical results are presented to evaluate the proposed LSA-SC system. In the 
simulation, we considered a quadrature-phase-shift-keying (QPSK) modulated SC signal with 
roll-off factor (3 = 0.25 and symbol duration T = 0.2 ps. The BS is equipped a ULA with 
M = 100 antennas without specification. A normalized extended typical urban (ETU) channel 
model [|T9ll is used in the simulation, which has 9 taps with maximum delay spread 5 ps. Both 
spatially independent and correlated channels are considered in the simulation. For practical 
implementation of the LSA-SC receiver, we adopt a twice oversampling based approach where 
the analog MF can be implemented in a discrete manner lilTTl . 

Fig. [5] shows SER versus SNR when channels are spatially correlated. From the figure, SER 
improves by increasing the array length since the tap gains corresponding to different antennas 
are less correlated for larger array length, which coincides with our analysis in Section IV. 

Fig. [6] shows SER versus array length. From the figure, SER can be improved by increasing 
the array length. When D = 50A, the SER can be hardly improved by further increasing the 
array length since the tap gains corresponding to different antennas are sufficiently independent. 
This corresponds to a distance of A/2 between adjacent antennas, which is exactly the Nyquist 
sample distance in space domain lf2TTl . 
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Fig. 7. SER versus antenna number with E s /Nq = —10 dB. 

Fig. |7j shows SER versus antenna numbers with E s /N 0 = —10 dB. From the figure, SER 
reduces as the increasing of the antenna numbers. However, severe performance degradation can 
be observed due to spatial correlation. 

VI. Conclusions 

We have investigated SC modulation for the LSA system in this paper. The LSA-SC receiver 
has been derived based on the waveform recovery theory, which can provide a straightforward 
way to understand the behavior of the SC receiver in the large-scale regime. When the antenna 
number is infinite and the channels at different antennas are independent, the desired waveform 
can be recovered except for the additive noise and the information-bearing symbols can be thus 
obtained by sampling the recovered waveform directly. We have also analyzed the non-ideal 
environments where the the antenna number is finite or the channels at different antennas are 
correlated. Our work shows that the LSA-SC is a good combination of LSA and SC, and it will 
be a competitive candidate for the uplink in future wireless systems. 
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VII. Derivation of E {\p(t )| 2 } 

Assuming the transmit symbols are independent, Pisi can be given by 

Pisi = £ E {|/[n]| 2 }. (A.l) 

n^O 

From (l29l) . E{|/[n]| 2 } can be expressed by 

E{I/MI 2 } 

= jp EEEE g(nT - r h + T pi )g(nT - r h + t P2 )- 

Pi h P2 h 
M M 

^^{ a mAh] a *m 1 \Pl} 0 ‘*m 2 h\ 0 ‘m 2 \P2}}- (A.2) 

17ll=l 7712 = 1 

Using the moment and cumulant relation Il22l . we have 

E {a mi [h]a* mi [pi]a* m2 [k\a m2 [pa]} = 

Cum{a mi [h], a* mi [p x ] , a* m2 [l 2 \ , a m2 [ p 2 ]} + E {a mi [h]a* mi [pi\}E{a* m2 [l 2 ]a m2 \p 2 }}+ 
E{a mi [h]a* m2 [l 2 ]}E{a* mi \pi]a m . 2 \p 2 }} + E {a mi [h]a m2 \p 2 ]}E{a* mi [pi}a* m2 [l 2 ]}. (A.3) 

Since a mi [l\], a* mi [pi], a* m2 [l 2 ], a m2 [p 2 \ are Gaussian distributed random variables, we have 

Cum{a mi [/i], a* mi [pi], a* 12 [l 2 \ , a m2 [p 2 ]} = 0, (A.4) 

and thus 

E{a mi [/i]a mi [pi]<r m2 p2]o m2 [P 2 ]} ~ 

E{a mi [h]a* mi \pi]}E{a* m2 [l 2 ]« m2 [p 2 ]}+ 

E{a mi [h]a* m2 [k)}E{a* mi [pi]a m2 [p 2 }}+ 

E{a mi [h\a m2 \p 2 ]}E{a* mi \pi]a* m2 [l 2 ]}. 


(A.5) 
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From (1241) . we have 

E{a mi [Zi]a^ 1 [pi]}E{Q!^ !l [Z2]a m2 [p2]} = **- Pi]<S[Z 2 - P2], 

E{a mi [Zi]oi^ 2 [Z 2 ]}E{Qi^ i [p 1 ]a' m2 [p 2 ]} = cx^p 2 (^JfZY ' y) “ Wpi ~ P2I 

Meanwhile, note that 


E{ct mi [/ 1 ]a m2 [p 2 ]} = E{«y [pijay [Z 2 ]} = 0. 

Substituting (IA. 6 I) to (IA. 8 I) into (IA.5I) . 

E {oimi [y a mi [Pl] tt rr!2 M Ctm2 [P 2 ]} 

=cr fi cr ia <5 [ Z i - pM 1 2 - P 2 ] + crla^p 2 ' t) ^ ~ ~ 

and thus 


E{l/N 2 } = 9 2 {nT) + tr ^ ^ a p a i g 2 (nT - r t + t p ), 

where we used the identity 


v 1 


M 


M 


EE ^ 2 

77li = 1 777-2 — 1 


m 1 — m 2 

M- 1 


D 

X 


= tr{R 2 }, 


for equation (IA.1QI) . 

Finally, substituting (IA. 101) into (IA. II) . Pisi is obtained as 

Pisi = t -^T~ <T p a i9 2 ( nT -n + T P ), 

n^O p l 

since g(nT) = 0 for n 7 ^ 0, which is exactly (l32l) . 
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